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■ Abstract 

0> I Heuristic insights into a physical picture of Davydov's soHtonic model of the 

• one-dimensional protein chain are presented supporting the idea of a 

^ ' non-equilibrium competition between the Davydov phase and a 

< : »„p,e™..„,. dy™„.cal- ■fa.oe.ee.nc- phase along the cha.n. 

■"nI" ' PACS 87.10.-)-e General, theoretical and mathematical biophysics 

(N ; 

. PACS 87.22.-q Physics of bioenergetic processes 

^ ' A well-known application of the cubic nonlinear Schrodinger equation 

(NLSE) in biology is Davydov's model of energy transport by means of soli- 
' tons along a-helix protein chains [0. Over the years, the model was turned 

. almost into a paradigm of bioenergetics |^ , although the experimental results 

^-H ' are debatable . 

, The a-helix is the most common secondary structure of proteins entailing 

' three spines in the longitudinal z-direction that we consider infinite in extent 

, and having the peptide sequence ( H — N — C = 0)„, n = oo, where the 

' dashed lines represent hydrogen bonds. There are side radicals and their order 

^ , is characteristic to each protein. The 3 spines are weaved together in a sort 

' of incommensurate quasi-one dimensional structure, but here we shall consider 

Ch , only the simple one-spine chain, containing the peptide groups as molecular 

units. 

In the Davydov model the main assumption is that about half of the chemical 
energy ( e ~ 0.422 eV) released in the hydrolysis of adenosine triphosphate 
(ATP) turns into the vibrational energy (eo = 0.205 eV) of the self-trapped 
' amide-I (C=0 stretching) mode of the peptide unit. The amplitude of the amide 

mode is self-trapped in the form of a sech- envelope soliton whenever there is 
a balance between the dipolar nearest-neighbour interaction and an admittedly 
quite strong nonlinear amide-phonon interaction. The energy transfer in the 
standard Davydov model is through the nonradiative resonant dipole- dipole 
interaction and a rather ambiguous vibrational- acoustic coherent state, the so- 
called Di ansatz. It is a mixture of quantum and classical Hamiltonian methods 
that has been deeply scrutinized in the literature Q . Davydov Hamiltonian can 
be written in the form 

Hd = Hc=o + Hph + Hint (1) 
Hc^o - J2'^^nBn-JiBiBn+i + Bl^,Bn) (2) 
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= J2hn,ib% + ^) (3) 

9 



N ■ 



^e-mail: rosu@ifug.ugto.mx 



1 



The capital and small operators are vibrational and phonon ones, respectively. 
Brown Q has shown in a clear way that Davydov Hamiltonian is a particular 
case of the general Frohlich Hamiltonian of polaronic systems |^ 

HfP = ^ Jfnnalnan + ^ flUJqb\bq + ^ hujq{xlb\ + xVbq)a\an (5) 
mn q qn 

with Di states satisfying the Schrodinger equation of the Frohlich Hamiltonian 
in the limit Jmn = and another D2 ansatz valid for Schrodinger evolution in 
the limit Xr — 0- 

The NLS subsonic soliton solutions of the energy transport coming out from 
the Davydov model are 



a{Cl = v/AV2e[^[^~^=*llcosh~i(-^C) (6) 

R 

PiO = Y-^^"'(^0 (7) 

m = ^(l-tanh(gO) (8) 

where ^ = x — Vst is the moving frame coordinate, J is the hopping (dipole- 
dipole) constant, x is the nonlinear dipole-phonon coupling parameter, 7^ = 
l/Vl — (s — Vs/va) is the soliton 'relativistic' factor, w is the elasticity 
constant of the chain, /i = x^ll/Jw, Pq = Q = MRx^ll/^wh^ , = 

eo — 2 J + h^v^/AJR^ — M is the peptide molecular mass, and R is the 

hydrogen-bond length. The first soliton is the vibrational one in which one may 
contemplate the soliton energy Eg self-trapped by the carrier wave. The second 
solution is related to the local deformation produced by the vibrational soliton 
in the lattice. The /3-kink is a domain-wall configuration of the displacements 
of the peptide groups into the new equilibrium positions. It has an enhanced 
stability of topological origin since all the peptide groups from the right side of 
the kink > ) are in nondisplaced positions whereas all the peptide groups at 

the left < 0) are displaced by the same amount (3o — order to destroy 

the domain wall configuration, one should first turn the left peptide groups to 
their initial position. Now let us state the essential point in our phenomenology 
which is to think of the slowly moving Davydov P-kink as an interphase boundary 
for the non- equilibrium transition from the Davydov dynamical regime of the 
polypeptide chain to a dynamic- 'ferroelectric' phase of the chain. This is not 
the only possible interpretation. Indeed, in the same 'non-relativistic' limit, the 
/3-kink, per se, is the solution of a classical anharmonic oscillator equation with 
a frictional term, but only for a particular value of the dissipation coefficient 
[0. However, here we shall pursue the first interpretation. In the literature 
on ferroelectricity it is common to consider the interfacial boundary as the 
kink solution of a time-dependent Ginzburg-Landau (GL) equation Q. Our 
interpretation is based on the fact that for low subsonic regime (s^ <C 1) the 
Davydov /3-kink is just the complement of the (quasi)-static double-well GL 
kink (i.e., Kjj cx (1 — Kql)). Since we are in a non-equilibrium situation 
the more precise terminology for these kinks is dynamic interphase boundaries 
or interfacial patterns. One would like to study their morphology during the 
growth. This is a difficult task since we are in a more complicated case as 
compared to the simple solid on solid model, or the kinetic Ising model, where 
the width of the mean field interface is given in terms of the nearest neighbor 
exchange interaction |^, corresponding to the J parameter in the Davydov 
model. The width W = 1/2Q of the /3-kink is determined by a number of 
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parameters, among which the J of nearest neighbors and the x parameter of the 
nonUnear interaction are controUng the interfacial morphology. It is also known 
that in the mean field theory the intrinsic width is identical to the correlation 
length in the direction perpendicular to the interface. 

Perhaps, one should also notice the formal analogy between the form of the 
Davydov kink and the Glauber transition rate in one-dimensional spin chains. 
The message of this analogy is that the kink is just a step structural function 
which is required by a Hamiltonian evolution and by a detailed balance condition 
in the spatial coordinate. 

The implicit occurence of the double-well GL kink in the Davydov /3 dis- 
placements offers the opportunity to introduce the multifractality issue and 
closely related intermittent features, which presumably are concepts of much 
interest for the dynamical phases of biological chains. Previously, Brax 
showed the equivalence of the GL equation with real coefhcients under random 
initial conditions and the linear heat equation with Gaussian random potential 
and provided a multifractal analysis of the problem of direct relevance for our 
considerations. As it is well-known the GL equation, which is a cubic reaction- 
diffusion equation, describes phenomenologically the evolution of the order pa- 
rameter in superconductive phase transitions, and models also spatial and time 
fluctuations of systems near Hopf (oscillatory) bifurcations. Moreover, in 1972, 
Scalapino, Sears and Ferrell ||l^ studied in detail the statistical mechanics of 
one-dimensional GL fields. They remarked that such fields can describe the 
dynamical behaviour of nearly-ordered systems which are not undergoing sharp 
phase transitions, and conjectured that the real-field case may have applica- 
tion in some organic chain systems. Our arguments support that idea. In the 
D'Alambert variable {£, = x — Vst), the GL kink Kcl is the solution of a GL 
equation of the type 



with real p and r coefficients. In the static case such GL equations are typical 
for the structural phase-transitions Q in equilibrium situations, but they can 
be used also in nonequilibrium, slowly-driven systems. Following Brax, see also 
[ p^ , one can develop a multifractal (MP) formalism, and also study the weak 
intermittent properties of the GL equation (and also of the Davydov's model), 
if the kink probabilistic distribution function is identified with the r- function 
of the MF formalism 



The g-parameter is well-known in the MF approach, being a convenient way of 
describing the mathematical properties of the sets of local growth probabilities of 
a physical quantity. We recall that in the MF formalism the function r(q) when 
it exists characterizes the power-law scaling of the cumulant generating function 
[^3|. The derivative of this function with respect to q is denoted by a and via 
the Lagrange multipliers procedure of statistical thermodynamics, one obtains 
the function f{a) which is the density of a measure, and can be interpreted as 
a fractal dimension when it is positive, and related to instabilities for negative 
values. The equations a = Sr/Sq and f{a) ^ qa — t represent the basis of the 
MF formalism, see e.g., |l^. Formally, q is the inverse temperature, r is the 
Gibbs free energy, and / is the entropy. One can plot the second derivative of 
the T-function (the 'specific heat') and find out intervals of the 'temperature' 
variable, within which the plot clearly displays features of a phase transition, 
that is a peak in the 'specific heat' at a certain value of the 'temperature' 
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variable p^ . At the same time, in our interpretation, q is the exponent of the 
local mean polarization of electric dipoles over sets of self-similar dipolar 
clusters, which in fact is another common assumption of the MF formalism. We 
hope to study in detail these heuristic claims in future publications. 

In conclusion, in our opinion, the old conjecture of Scalapino, Sears and 
Ferrell concerning incipient GL-type phase transitions in one-dimensional 
organic chains might be at a quite good place within the dynamical Davydov's 
model of protein chains. Moreover, one can attribute self-organized features to 
the non-equilibrium competition between the two dynamical phases as Canessa 
and myself have shown in a previous paper IT^ . Multifractal and intermittent 
scaling in the present biological context is of great interest. 

Regarding our conjecture of relating Davydov's model to dynamic-structural 
transitions, it is to be also mentioned the debate on the connection between the 
equations of motion for Davydov's solitons and the (/)^-chain [^ . 

Finally, we recall that metastable, biological- 'ferroelectric' states have been 
first considered by Frohlich 
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